arXiv:cond-mat/0311604vl [cond-mat.str-el] 26 Nov 2003 


A renormalized Hamiltonian approach to a resonant valence bond wavefunction 1 

F C Zhang, C Gros, T M Rice and H Shiba 2 
Theoretische Physik, ETH-Honggerberg, CH 8093 Zurich, Switzerland 

The effective Hamiltonian of strongly correlated electrons on a square lattice is replaced by a renor- 
malised Hamiltonian and the factors that renormalize the kinetic energy of holes and the Heisenberg 
spin-spin coupling are calculated using a Gutzwiller approximation scheme. The accuracy of this 
renormalization procedure is tested numerically and found to be qualitatively excellent. Within the 
scheme a resonant valence bond (RVB) wavefunction is found at half-filling to be lower in energy 
than the antiferromagnetic state. If the wavefunction is expressed in fermion operators, local SU(2) 
and U(l) invariance leads to a redundancy in the representation. The introduction of holes removes 
these local invariances and we find that a d-wave RVB state is lowest in energy. This state has a 
superconducting order parameter whose amplitude is linear in the density of holes. 


I. INTRODUCTION 

Very soon after the discovery of high-T c , superconduc- 
tivity Anderson [1] proposed that it was caused by a co- 
operative condensation of carriers moving in a resonant 
valence bond (RVB) state of spins. Since then this pro- 
posal has been studied extensively and the best account 
is in Anderson’s recent lecture notes [2]. Many other pro- 
posals have been made (for a Review see [3]), and there 
have been questions raised about the use of the simplified 
effective Hamiltonian derived from a single-band Hub- 
bard model in the atomic limit that forms the starting 
point of Anderson’s treatment. We shall not go into these 
questions here but just point out that two of us (Zhang 
and Rice [4]) recently gave an explicit demonstration that 
a two-band model describing hybridised copper 3d and 
oxygen 2p states can also be reduced to the same effective 
Hamiltonian in an appropriate limit. 

The effective Hamiltonian contains the strict local con- 
straint which forbids double occupancy of any site. This 
constraint is very difficult to handle analytically. One 
of the most physically transparent methods to treat this 
type of problem analytically has been the Gutzwiller ap- 
proximation which introduces a renormalization of the 
quantum mechanical expectation values by a classical 
weighting factor [5]. Such renormalisation can then be 
incorporated into a Hamiltonian which may be treated 
by conventional methods. This approach was used for 
the heavy fermion problem by Rice and Ueda [6] and 
was shown to be equivalent to an optimal slave-boson for- 
mulation by Kotliar and Ruckenstein [7]. In this paper 
wc will consider this renormilisation Hamiltonian method 
for the effective Hamiltonian. Unlike the renormalised 
Anderson Hamiltonian studied by Rice and Ueda [6], in 
the present case the renornralised Hamiltonian cannot 
be simply diagonalised and we must resort to a further 
mean-field approximation. Mean-field approaches have 
been considered by many authors [8], [9], [10]. Here we 
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choose to formulate the problem in terms of a varia- 
tional wave- function. This has several advantages. First 
it shows us that a consistent mean-field theory must be 
formulated in terms of two expectation values i.e. one 
must include particle-hole amplitudes of the form (c + c) 
in addition to particle-particle amplitudes of the form 
(c+c+). This point has been recently realised by others 
as well [11] + . The coupled equations to minimise the 
energy have a wide class of degenerate solutions at half- 
filling. Secondly a wave-function formulation is suited to 
examining the role of the redundancy in the fermion rep- 
resentation which is not present in the spin representa- 
tion. At half-filling, this redundancy which has its origin 
in the reduction from 4 to 2 degrees of freedom per site 
as one goes from fermion to spin representation, appears 
as a local particle-hole (517(2)) and gauge (U( 1)) invari- 
ance. The large degeneracy of the mean-field description 
arises from this redundancy and it can be shown that it 
corresponds to the same state in the spin representation. 
Further the appearance of coherence in the fermion rep- 
resentation is illusory so that there can be no true phase 
coherence as stressed by Baskaran and Anderson [12], 
Thirdly, this formulation allows a direct comparison with 
the variational Monte Carlo (VMC) results. This allows 
us on the one hand to check the validity of the renor- 
malised mean-field theory and on the other hand it gives 
us more insight into the numerical VMC results. Both 
qualitatively and even quantitatively good agreement is 
found; for example both point to a d-wave paired state as 
the most stable and a true superconductivity order pa- 
rameter which vanishes at half-filling and grows linearly 
in the deviation from half-filling. The largest discrepancy 
occurs for the antiferromagnetic state which, within this 
scheme is higher in energy than the RVB state, contrary 
to the VMC results. 

Our treatment is essentially limited to zero tempera- 
ture and the extension to finite temperature will be non- 
trivial. Some discussion of the problems of calculating 
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excitation energies is given. In particular there are two 
energy scales of excitations given by the gauge coherence 
energy (determined by the kinetic energy) and the mag- 
netic coherence energy respectively. Anderson[2] has em- 
phasised this splitting of the charged excitations (holons) 
and spin excitations (spinons). 

II. THE MODEL AND THE RENORMALISED 
HAMLLTONLAN 

We study the Hubbard model on a square lattice. In 
the limit of large on-site Coulomb repulsion U and at one- 
half, or, slightly less, filling the Hubbard Hamiltonian can 
be transformed to the form 

H = H t + H s (1) 

Ht = —t ^2 C t,<T C j,V + HC 

H s = 

(hj) 

with the local constraint the number of electrons on any 
site <1. This transformation has a long history, and has 
been used by [13], amongst others. In (1) H t , and H s , 
are the kinetic and magnetic energies respectively and 
(i,j) represent the nearest-neighbour pairs. S) are the 
spin = | operators and J = . We neglect terms which 

are higher order in the small parameters t/U and the 
hole concentration 5 (= 1 — n where n is the electron 
concentration) . 

Since the high-T c , superconducting materials show 
strong antiferromagnetic (AF) spin correlations [14] it 
is believed that this model contains the essential physics 
for the high-T),, superconductivity [1] , [4] . 

To study the ground state and the excited states of (1), 
we use a projected BCS trial wave-function as suggested 
by Anderson [1] for a RVB state: 

\V>) = p d l^o) ( 2 ) 

IWo) = IT( Ufe + VkC k'\ C -kO 1°) ( 3 ) 

k 

where the Gut z wilier projection operator Pd = na- 

i 

n^riii) and |0) is the vacuum state. Uk and Vk are the 
variational parameters satisfying the normalization con- 
dition for |y> 0 ): \uk\ + \vk\ = 1. 

Some special forms of (2) have recently been studied 
numerically. Using the VMC technique, which treats the 
projection operator exactly [15], [16], [17], [18], [19], the 
energies of these states have been numerically calculated. 
It is found that in the square lattice, the projected Fermi 
liquid state (i.e. the state with u k v* k = 0) is unstable 
against d-wave pairing [16]. At half-filling the energy of 
the d-wave state is found [17] to be very close to the 
ground-state energy extrapolated from the exact small 


system calculations [20]. In contrast to the extrapolated 
exact small system calculations, the d-wave trial wave- 
function has no long range antiferromagnetic order [17] 
and may therefore be viewed as an example of a quantum 
spin liquid. A VMC study of superconductivity has been 
made independently by Yokoyama and Shiba [19]. They 
also concluded a possibility of a d-wave superconductiv- 
ity. 

The projected BCS wave-function is a natural gener- 
alization of the usual BCS state to strongly correlated 
systems. The projection operator, however, makes diffi- 
culties for an analytic approach. In this paper we shall 
use a renormalised Hamiltonian approach to treat the 
projection operator, and systematically investigate the 
state (2), carrying out explicitly the variational proce- 
dure. In this approach following Gutzwiller [5] the effect 
of the projection operator on the doubly occupied sites 
is taken into account by a classical statistical weighting 
factor which multiplies the quantum coherent result cal- 
culated with |y> 0 ). A clear description of the method has 
been given by Vollhardt [21]. The hopping energy and 
the spin-spin correlation of the nearest neighbour sites in 
the state | ip) are related to those in the state | <p 0 ) by 

( c tr c 3<r) = 9t(c+C ja ) 0 

(Si ■ Sj) = gs (Si ■ Sj) Q 

where (A) 0 , and (A) are the expectation values in the 
states |^o), and \<p) respectively. The renormalization 
factors gt and gs are determined by the ratios of the 
probabilities of the corresponding physical processes in 
the states \<p) and |y> 0 ). In figure 1 we illustrate the pos- 
sible hopping processes in these two states. The proba- 
bility of such a process in the state | <p) is 

[njT( 1 -n i )n iT ( 1 -n j )] 1/2 
while that in the state |<p 0 ) is 

[r 1 j T ( 1 — n d) n *T ( 1 — n il )] 1/2 

n-ia- are the average electron occupation numbers, (n, = 
Yh n iv) which are the same in the states | ip) and |<p 0 ), 

a 

because of the spin symmetry of the wave- functions. This 
leads to the result [21] 

gt. = 25/ (1 + S) (4a) 

To determine gs, we consider the spin exchange process 
shown in figure 2. The spin-spin interaction occurs only 
when both sites are singly occupied. The probability for 
such a process in the state \tp) is (njfnjj,n J qnjT) 1 / 2 , while 
in the state |y> 0 ) it is [njf(l — 1 — nq)n,q( 1 — 

n jt) n d(l — n a )] 1/2 . The same result is obtained for the 
z component interaction S?Sj. Thus one finds + 

g s =A/(l + 6) 2 (4b) 


[4] +In a systematic series expansion on S and J/t , the higher-order 
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FIG. 1: The possible hopping processes (a) in the non- 
projected pairing state (3) and (b) in the projected BCS state 
(2). The spins with broken arrows are optional in the (a) 
configurations. 
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FIG. 2: The spin exchange process in the states (2) and (3) 

It is important to realise that the projection operator 
greatly enhances the spin-spin correlations. To further 
illustrate this point, we list in table 1 all the possible 
two-site states together with their weights and the con- 
tributions to the spin-spin correlation in the half-filled 
case. 

Having determined the renormalisation factors, we can 
define a renormalised Hamiltonian given by 

H'=g t H t +g s H s (5) 

The energy of the system in the state \<p) can be evaluated 
as the expectation value of H in the state |<p 0 ) 

W = ( H ‘) o < 61 

Equations (4a,4b)-(6) form the basis of our renornralised 
Hamiltonian approach, which is analogous to the ap- 
proach used by Rice and Ueda [6] for the periodic An- 
derson model with the difference that here we make a 
further mean-field approximation. This is because in the 
context of the periodic Anderson Hamiltonian the most 


terms of S in (4a, 4b) should be dropped away to be consistent 
with the effective Hamiltonian (1), where the higher-order terms 
are not included. This, however, does not change the qualitative 
physics discussed in this paper. 
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FIG. 3: This table illustrates the enhancement of the spin-spin 
Correlation in the projected BSC state \ip), equation (2), over 
that in the BCS state |<p} 0 (3), at half-filling. The weight of 
the configurations actually contributing to (Si ■ Sj) increases 
by a factor of four due to the projection. The configurations at 
each site are denoted by 0 (empty state), Tj (doubly occupied 
state), and a (singly occupied state with spin a) 


important physical effect is the renormalisation of the f- 
level to the Fermi surface and not the spin interaction, 
which would make an exact treatment of the effective 
Hamiltonian impossible. 

To justify this approach, we have carried out Monte 
Carlo calculations. Figures 3-5 show the comparisons 
between the renormalised mean-field theory and the es- 
sentially exact MC results for these wave- functions. The 
quantitative agreement is within 5-15%, while the quali- 
tative agreement is excellent for the wave- function (2). 

After replacing the projection operator, the energy of 
the system can be evaluated analytically. The variational 
task is to minimise W in (6). This leads to coupled gap 
equations, which we will derive and solve in the following 
sections. 


III. GAP EQUATION 

In this section we derive the gap equations for the 
projected BCS wave-function within the renornralised 
Hamiltonian scheme described in section2. We consider 
only the even-parity case, i.e. , U-kV*_ k = uyv), and 

\v-kf = \v k \ 2 . 

Evaluating (6), we obtain 
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FIG. 4: A comparison between the renormalised mean-field 
theory (RMF) (see fa, fb, 5, 6) and the Monte Carlo (MC) 
result for the kinetic energy ( T ) per hole in the projected BCS 
state (2). Both were calculated with a total number of sites 
Ns = 82 and a number of holes Nh = 8, 16. The varia- 
tional parameter A is related to the parameters of the state 
(2) by tit. — ' 0,1 /■> where u n is a param- 

eter, and Ek is given by (8b). In the d-wave pairing state, 
A*, = A(cos(k x ) — cos(ky)) and in the s-wave state, Ak = A. 
The full circles and squares are the MC results for the s- and 
d-waves respectively. The dotted and broken curves through 
the MC results are guides for the eyes. The second pair of 
dotted and broken curves are the results from RMF for s- and 
d-waves respectively. 


W = 2 Strife M 2 (7) 

k 

+^E(WM 2 \ V k '\ 2 + U kV* k V k -U* k ,) 

k,k' 

where Ns is the total number of sites, and 

3 

Vfc = --gsJlk (8a) 



FIG. 5: A comparison between the renormalised mean-field 
theory (fa, fb, 5, 6) (RMF) and Monte Carlo (MC) results 
for the nearest neighbour spin-spin correlation (Si ■ Sj) in the 
projected BCS state (2). Both were calculated with a total 
number of sites Ns = 82 and a number of holes Nh = 0, 8, 16. 
The variational parameter A is related to the parameters at 
the state (2) by where g 0 is a 

parameter and Ek is given by (M). In the d-wave pairing state, 
Ak = A(cos (k x ) — co s(k y )) and in the s-wave state A k = A. 
The full circles and squares are the MC results for the s- and 
d-waves respectively. The dotted and broken curves through 
them are guides for the eyes. The second pair of dotted and 
broken curves are the results from RMF for the s- and d-waves 
respectively. 


Let pi be the chemical potential of the system, the 
quantity we want to minimise is 

W = (h' - piN^ 


= -tlk 

(8b) 

7 k = 2 (cos k x + cos k y ) 

(8c) 


Note that Ek and Vk have the same functional form, 
since H$ is derived by kinetic exchange. The electron 
number operator N = Yl c k a Ck ^ -^ as expectation value 

k,a 

<7V) = 2£M 2 . 

k 


While minimising W with respect to Uk and Vp for fixed 
H, one must remember that H is also a function of 5 due 
to the renormalisation factors, and hence a function also 
of Vk- Carrying out this procedure, we find that 


Kl 2 

kl 2 


-[i + — ] 

2 E k J 

2 E k J 


(9a) 


(9b) 
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FIG. 6: The nearest neighbour spin-spin correlation function 
(Si ■ Sj) as a function of electron filling in the projected Fermi 
liquid state (ukVk = 0 in (2)). The renormalised mean-field 
theory (ja, jb, 5, 6) (RMF), broken curve, agrees well with 
the Monte Carlo (MC), full circles, result in the entire filling 
region. 


with t = x and y, i + r denotes the NN of i in the r 
direction. Since we consider the even-parity case, C T is 
real, but A T can be complex. A T and £ r , satisfy the 
following coupled gap equations: 

Afc = N s 1 E T k-k' 

k’ k 

C k =e k + Ns 1 Y. lk _ k ,^- ( 14 ) 

k 

The first one is the same as the usual BCS gap equation. 
The second one originates from the particle-hole corre- 
lation. From (12), it is clear that A k is related to the 
pairing in the unprojected state |<^ 0 ). It describes the 
’smearing’ of the pseudo-Fermi surface. However, A*, 
is not the superconducting order parameter in the pro- 
jected state | ip) in our theory. E k turns out to be the 
quasi-particle excitation energy (in units of jgsJ) in the 
pairing state as we will show in section5. 

The coupled gap equations (14) are the basic equations 
in our approach. They can also be written in the x and 
y component form: 


and 


u k vl = A fe /2 E k (9c) 


= N s cos hr 

S ^E k 


where 

E k = (Cl + A 2) 1 / 2 (9d) 

The parameters A^and C k are dimensionless, and they 
are related to the particle-particle and particle-hole pair- 
ing amplitudes respectively 

A fc = Aa, cos (k x ) + A y cos (k y ) 

(k=£k~ Cx cos (k x ) - C y cos (k y ) (10) 

where 


Cr = -^S 1 E^COSfc T (15) 

k 

The gap equations must be solved simultaneously with 
the hole concentration equation, 6 = 

k 

Before we discuss the non-trivial solutions, we note 
that Afc = 0 is a trivial solution of the gap equations. 
This corresponds to the projected Fermi-liquid state. In 
this case, C k changes sign at the Fermi surface. The pa- 
rameters C x = C y { = C) are given by 

C = Ng 1 ^2 cos k x + cos k y 

Cfc<o 


~ _ ( 9t£k ~ k) 

£k — 3 T 

jdsJ 

and Ji is related to the chemical potential y by 


The volume of the Fermi sea is determined by the number 
of electrons. The energy per site is 

w = -4, gt tc- ^ gsJ( 2 


In (10) 



In particular, w = — 48/7T 4 J ~ — 0.49J in the half-filled 
(11) case. It will be shown in the next section that this trivial 
solution is unstable against the pairing states with Afc ^ 

0. 


A t 


,+ „+ 

'iT C i+rl 


- c u c trr) c 


(12) 


IV. SOLUTIONS TO THE GAP 
EQUATIONS — HALF-FILLED CASE 


Cr — E ( C ta C i+rcr) 0 
a 


( 13 ) 


At the half-filling , S = 0, y = 0 and there is no kinetic 
energy. We are interested in the possible lowest energy 




states of the solution of (14). The total energy of the 
system has a simple form in this case by use of (1) and 
(14) 
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w = --gsJ^Ek 
k 

Therefore the lowest energy states correspond to the 
maximum value of ^Ek- For this reason, we use 
k 

an ansafz for Ek to examine the solutions of the gap 
equations+ 


E k = C(cos 2 (k x ) + cos 2 (fc y )) 1/2 (16) 

where C is a parameter to be determined. Note that 
such a choice gives only four point zeros for Ek . The gap 
equations then reduce to a single equation, and we get 

C = ^^( cos2 ( fca 0 + cos 2 {k v )) 1/2 

which has numerical value C « 0.479. The energy per 
site is 

w = ~ g s JC 2 = -0.688J 

This energy is much lower (about 20%) than that of the 
projected Fermi-liquid state found above. The parameter 
A k describes the pairing correlation in the renormalised 
Hamiltonian. Finite values of A k indicate the binding of 
the electron pairs in the pairing states (3). 

We now determine the parameters A T and C T , required 
for the choice of Ek in (16). Using (9d) and (10), we 
find that they should satisfy the following simultaneous 
equations: 


Cx + 




= c 2 


2( x ( y + (A x A* y + HC) = 0 (17) 

There is a wide class of parameters which satisfy the 
conditions (17). All the states in this class give the same 
expectation value of the renormalised Hamiltonian H in 
(5). Therefore at the half- filling, H has a large degen- 
eracy of ground states. For real A^ and A y . these states 
can be illustrated diagrammatically as shown in figure 
6(a). A few examples of these states are 
d-wave pairing : 

A* = -Ay = ( x = ( y = C/V2 (18a) 


[5] 4 .There might be other solutions for the gap equations. We be- 
lieve that theform of the energy given by (16) givesthe lowest 
energy. 



FIG. 7: (a) Diagrammatic illustration of the degenerate 

ground states for the renormalised Hamiltonian (5) at half- 
filling. C, = + C y V> A = A x x + A y y with ( T , A T given 

by (12) and (13). The full arrows represent the d-wave pair- 
ing state with ( i A, and |C| = |A|. All the states in (17) 
with the real parameters A T may be obtained by rotating ( 
and A simultaneously by angles 6, or reflecting the two vec- 
tors about the x axis, (b) An illustration of the SU( 2) gauge 
transformation. A state described by (( ,A ) in (a) can be 
obtained by a local SU( 2) from the state ((, A), under which 
cf„ at the four sites of the plaquette transform according to 
(19), with (onfli) as denoted. The minus sign in front of the 
parentheses corresponds to the states in (a) after a reflection 
about the x axis. The transformation operator a tCT at other 
lattice sites is determined by a translation. 


d-wave density matrix: 

A x = A y = C x = -Cy = C/V 2 (18b) 

chiral state: 

A* = -iA y = C, C* = C y = 0 (18c) 

anisotropic state: 

A, = Cy = C, Ay = C x = 0 (18d) 

We remark that the d-wave density matrix state is 
different from the extended s-wave state proposed by 
Baskaran, Zou and Anderson [8]. In their theory, the 
particle-hole amplitude (c + c) 0 is not included, i.e. ( x = 
C y = 0. Therefore their state has the same energy as the 
projected Fermi-liquid state. The d-wave pairing state 
was studied numerically in [17] and [18], and the chiral 
state was discussed in [10]. They belong to the solutions 
of the same gap equations in the present approach. 

The degeneracy of the ground states of H may be 
explained using the local SU (2) symmetry of the Heisen- 
berg Hamiltonian Hg- This symmetry has been pointed 
out by Anderson [2] . Very recently Affeck and co-workers 
[22] have studied the invariance under a time-dependent 
SU( 2) gauge transformation, and discussed the equiva- 
lence of some different mean- field theories. Here we wish 
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to show that the symmetry described in (17) is a sub- 
group of the local SU(2), which preserves translational 
invariance and even parity that we imposed in deriving 
the gap equations. 

Consider a local SU(2) gauge transformation, under 
which electron operators at site i transform as 


At half-filling, each site is singly occupied. Therefore 
any half- filled state | if) transforms into itself except for 
an overall phase factor under the 517 ( 2) operator U: 

\<p)^U\v) = e ie \v),0 = Y J 0i 


E T 


+ Pi°il 


Cil ->■ + OL*Cii (19) 

where oaa* + /3,/3* = 1. These are the particle- hole 
transformations with spin conservation. Hg, hence H 
at half-filling is invariant under these transformations. 
Therefore all the states related by (19) are degenerate. 
A set of the parameters A T and </ T , in our gap equa- 
tions transforms to another set of the parameters under 
(19). The transformations corresponding to the degen- 
erate state (17) derived from the gap equations are rep- 
resented in figure 6(b), for real values of A r . There is 
also one to one correspondence to transformations of the 
d-wave pairing state to a state with complex A T such as 
the chiral state (18c) with a and (3 complex. 

Local 17(1) gauge symmetry is a sub-group of the local 
SU (2) symmetry since it is of the form 



A general choice of 6 t transforms the BCS pairing state 
(3) to a non-BCS-type state, which has the same energy 
as the former. Since under such a transformation the 
Bloch coherence is lost, i.e. a state of the form 


Ck,( 7 


\ ' pik-Rj ~ . _ \ ' ik-Rj+idi 

3 3 


is no longer a coherent superposition of the original 
states. Yet we can equally well pair with c/ k , , 
and the energy would be the same as if we pair c k , with 
c! fc , . Thus the apparent coherent k-space pairing in 
the BCS wave-function (3) is illusory. The absence of a 
coherent pairing order parameter at half-filling as a con- 
sequence of the 17(1) gauge invariance has been stressed 
by Baskaran and Anderson [12]. 

It is important to realise that the states that are degen- 
erate due to the SU (2) gauge symmetry are the unpro- 
jected states |<p 0 ) of (3), rather than the physical states 
| (p) which obey the strict local constraint. How does a 
projected state change under 517(2)? According to (19), 
a vacuum state (empty state) at site i must transform 
under SU (2) as 

|0) i -e^(a^/3:c+c+)|0) i 


This ensures that the vanishing of the state c, )<T |0) i re- 
mains unchanged under the transformation as required 
physically. However a singly occupied electron state 
transforms under 517(2) as 


10), 


e i6i c + 

C S,<7 


|0>, 


Although U does not commute with the projection op- 
erator Pd, we observe for the half-filled state |y> 0 ) 

UPdW) 0 = PdU\y) 0 

One way to see this is to notice that there is no empty 
site state also in Pd |<p) 0 . Thus we can rewrite 

Pd W ) 0 = IIKt - n u ) 2 IWo 


The SU (2) transformations all commute with the opera- 
tor (n, iT - n,q) 2 . Let | ip) = P d |y>) 0 , and 
then 


<p ; o = u\<p) 0 , 



P d U |<p} 0 = U \ip) 


e ie W) 


e iB Pd \<p) 0 


This proves explicitly that the two states |y>) 0 and 

related by any local 517(2) gauge transformation corre- 
spond to the same state | ip) except for a phase factor. 
Therefore all the states in (17) correspond to the same 
physical state. The RVB ground state is non-degradable. 

There are also redundancies in the higher energy states 
in the fermion representation. For instance, the state of 
Baskaran, Zou and Anderson [8] at half-filled is identical 
to the projected Fermi liquid state, because the former 
transforms to the latter under (19) with a = 1,(3 = 0 
in one sublattice, and a = 0, /3 = 1 in the other. This 
equivalence was also pointed out by Yokoyama and Shiba 
[18] in a different way. 

We now comment on the local gauge symmetry in the 
original Hubbard model, which in terms of the original 
fermion operator is 



H h = —t E <<a- + HC + ' E^'E'''-. 

(*>.?) ic i 


This Hamiltonian is not invariant under local gauge 
transformations with respect to the operators d, jCr . How- 
ever, up to any finite order in t/U there exists a canoni- 
cal transformation, which eliminates the doubly occupied 
sites [13] 

H h -> e iS H H e~ iS 

At half-filling, all the odd order terms in t/U vanish. 
The Hamiltonian is locally gauge invariant with respect 
to the electron operators in the new representation, i.e. 
the c, jCT of (1) are Wannier operators of the old represen- 
tation, dj jCT . Therefore the local gauge symmetry holds 



to any finite order in perturbation theory in t/U. This 
is the same as saying that the system has undergone a 
transition to a Mott insulator [13]. 

We have so far only examined the projected BCS-type 
trial wave- functions. It is likely that the true ground 
state of the model Hamiltonian (1) at half- filled is the AF 
state [19]. Recently, Yokoyama and Shiba [18], [19] have 
studied a projected Hartree-Fock-type AF state. Using 
VMC they found the energy per site at half-filled to be 
— 0.642 J, slightly lower than — 0.636J, the value found by 
Gros [17] in the d-wave pairing state by using a similar 
technique. However we may argue that the holes favour 
the pairing state away from half-filled because of the gain 
in kinetic energy. We have also applied the renormalised 
Hamiltonian approach to the AF states. Within this ap- 
proximation, we find that at the half-filled, the AF state 
has higher energy than the RVB state (2), in contrast 
with the VMC results. We present the derivation and 
the results in Appendix 1. 


V. NON-HALF-FILLED CASE 
A. Ground state 

Firstly we examine the energy needed to create prop- 
agating Bloch states. The simplest states for holes have 
the form 


|^i,cr) — G,cr | 

which destroys a real electron at site i. We may also 
make a propagating Bloch state for the hole of the form 


l^p.cr) — Cp,(j | T 3 ) (20) 


to remove an electron and create such a Bloch hole is 
— 0.32t — 4aJ. 

The Bloch states are not the lowest energy states of the 
holes however. We now apply the gap equations to study 
a system with a few pair of holes. The energy to create a 
pair of holes is — 2/x by the definition of the chemical po- 
tential. Since the parameter /I = 0 at the half- filled, (11) 

gives the energy per hole to be Ng 1 , a quantity 

related to the unprojected state |<p) 0 at half-filling. In 
the presence of holes, the kinetic part of the Hamiltonian 
explicitly breaks the SU( 2) and U(l) gauge symmetries, 
while the Heisenberg spin part remains invariant under 
these symmetries. 

Using (4a, 4b) and (5), for the states |<p) 0 described by 
(17), the magnetic energy per hole is 


w s w < H *>° = 6C2J 


a value equivalent to the loss of four bonds in the spinspin 
correlations, and it is the same for all these states as a 
consequence of the SU(2) gauge invariance of the spin 
part of the Hamiltonian Hg. The kinetic energy per hole 
in this case is given by 

C T , is the particle-hole correlation in |<p ) 0 as defined in 
(13). When the holes are introduced, a fraction (g t ) of 
this correlation becomes coherent in the state | ip). There- 
fore the larger values of (j T , correspond to the lower ki- 
netic energy of the holes. But ( T , are subject to (17). 
The kinetic energy can be written in the following form 
by using (17) : 


T = —At(2C 2 


A X + A y 


) 1/2 sgn((( x +( y )) 


A rigorous calculation is possible at half-filling. 

We consider any translationally invariant spin singlet 
state \ip) at half- filling. Let us denote by a the NN spin- 
spin correlation in | ip), a = (Si ■ Sj). Then the magnetic 
energy loss of a hole in the state is — 4aJ, because 
the four bonds connecting the site i are mixing. Since 
the matrix of Hg in (2) between any states |$j i(T ) and 
|$j i<7 ) is diagonal, the moving hole state of (20) has the 
same magnetic energy as in | < f>i i(T ). 

The kinetic energy of the hole in (20) is given by 

(H t ) p = 2 (ni^nj t<T + <S% S~ exp(z p ■ R oi ) + HC 

(id) 

where ( ) denotes the expectation value in the half- 

filled state. Using the fact that (■ riirij ) = 1, we get 

(H t ) p = t( 1 + 4a)(cos(p x ) + cos (p y )) (21) 

Since a ~ —0.33 for the ground state, (21) gives a band 
width for a Bloch hole of 0.64 |f|. The minimum energy 


Different parameters A r , and ( T , describe different states 
with different energies upon doping. The above hole 
kinetic energy expression immediately leads to the im- 
portant conclusion that the d-wave pairing state, where 
A x + A y = 0 gives the best kinetic energy, which is 

T = -\\f2Ct = -2.71 1 

Both the d-wave density matrix state (18b) and the chi- 
ral state (18c) have zero kinetic energies, and are not 
favoured upon doping. The kinetic energy for the d-wave 
pairing state in our analytical approach is quite close to 
the VMC result, where it is found to be —2.55 1 for sys- 
tems with 10% holes [17]. It is also substantially below 
the value found for a Bloch hole. 

The introduction of some holes breaks the local gauge 
symmetry and causes the ground state to be coherent. 
The stable lowest energy state upon doping is the d-wave 
pairing state. 

It is worthwhile remarking that in a Hubbard model it 
is the expectation value of the kinetic energy which deter- 
mines the integrated optical weight associated with the 
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FIG. 8: Variational parameter A and superconducting order 
parameter A sc as functions of the hole concentration S for a 
choice of t/J = 5 in the d-wave pairing state. 

charge carriers in the f-surn rule [23] and in the present 
case this optical weight is proportional to the number of 
holes with an optical mass determined from the propor- 
tionality constant of order at -1 (where a is the lattice 
constant). 

The gap equations for the finite hole concentrations 
can be solved numerically. Here we shall consider only 
the most stable state-the d-wave pairing state. In this 
case, we set = — A y = A, and £ x = £ = £. The four 
equations (13) reduce to two because of the symmetry 
between the x and y components. These two equations 
uniquely determine A and £ for the fixed values of t, 
J and Ji. The numerical results of the gap equations 
are plotted in figure 7 for A as a function of the hole 
concentration. 

We now discuss the superconducting order parameter. 
As mentioned in section3, A is not the order parameter. 
The superconducting order parameter is 

A sc(Rij) = (<T c hi - c Mt) 

an expectation value in the projected state (2). This 
quantity describes the Cooper pairing in a real space rep- 
resentation. We shall adopt the Gutzwiller method to 
calculate this quantity. In analogy to the derivation for 
the hopping energy in section2 we find that the nearest- 
neighbour sites i and j 

( c tAi) = 9t ( c tAi) 0 

Therefore for nearest-neighbour sites, the order parame- 
ter is related to the variational parameter a in the gap 
equations by 

A sc = 9t A (22) 


The value of A sc as a function of S is plotted in figure 7 
in comparison with A. A sc vanishes linearly near (5 = 0. 
A sc found in our theory is in good agreement with the 
Monte Carlo results [[17]]. The absence of the supercon- 
ducting order parameter at half-filled obtained from (22) 
agrees with the discussion in section4 from the viewpoint 
of the local gauge symmetry. 

The kinetic energy of holes in the AF state is found to 
be quite high in our analytic approach (see Appendix 1). 
However, the Gutzwiller approximation we adopted is too 
rough to determine whether a RVB or AF state has lower 
energy. Numerical results of VMC [17], [18], [19] suggest 
both the spin-spin correlation energy and the kinetic en- 
ergy of the holes between the d-wave pairing state and 
the AF states are very close. The question which state is 
more favourable in energy remains unresolved. 

B. Excited states and finite temperatures 

We begin by examining the spin degrees of freedom 
in half-filled and near half-filled cases. An excited state 
can be created by applying the spin raising operator to 
a specific site to obtain 

\^ it+ ) = SfP d \<p) 0 

We can commute Sf with Pd, to obtain 
i*i,+> = PdSt \<p ) 0 

= Y. ei(p ~ p)Rzp AA c P'A) 0 

p,p 

This state is therefore a superposition of two independent 
quasi-particle states similar to a metal where the low en- 
ergy excitations are made up of superpositions of electron 
and hole states. The quasi-particle states (spinons) can 
be defined by 

I’hpd) = i°) ( 23 ) 

k^p 

The quasi-particle energy E p , is defined to be the differ- 
ence of the expectation values of A' = H — gN in the 
state Tpj) and in the ground state | ip). We use the 
Gutzwiller method to calculate the energy of the state 
(23). The energy difference between the two states con- 
tains two parts. One is due to the changes of the renor- 
malization factors gt. and g s , the other comes from the 
change of the wave-function itself. The former just can- 
cels exactly the second term in fj in (11). Using (1) to 
calculate the energy difference due to the wave-function 
change, we get 

= (l-2Vp){gt£ P + Ng 1 J2 V k -p v k-T i ') 

k 

-\-2u p v p Ng _1 E V k -pU* k v k 

k 
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Applying the gap equations to simplify the expression, 
we obtain 


3 

= —gsJEp 


(24) 


Note this energy is independent of the local SU( 2) gauge 
and does not depend on the particular fermion represen- 
tation. At the pseudo-Fermi surface, where by definition 


C = 0, we have E p 


\ 


. Since the state (23) breaks a 


pair of electrons, E p describes the binding energy of the 
pair at the pseudo-Fermi surface. The excitation energy 
depends on the particular RVB. For the Fermi liquid state 
(ukVk = 0) then E p = 0 over the whole pseudo-Fermi 
surface. However in the ground RVB state it vanishes 
only at four points, e.g. when n = 1, E p = E p = 0, if 
( PxiPy ) = (±7t/2, ±tt/ 2) and the density of spinon states 
at low energies is 


N{E) = Y, 5 ( E ~ E v) 


2 E 


<\9sJ) 2 ' 


as E — > 0 


We turn now to a brief discussion of the system at finite 
temperature. The extension of the mean-field gap equa- 
tions to finite T is not so straightforward. The existence 
of a finite A*, is controlled by the energy scale of Ek i.e. 
by J in (24). On the other hand if we consider the limit 
5 <C 1 there is a very small energy scale oc St which con- 
trols the definition of a coherent gauge. In other words 
it is only the kinetic energy which allows us to determine 
the gauge uniquely and at temperatures J T St, 
the gauge coherence will be lost. Yet in this temperature 
range the magnetic coherence survives since as we have 
stressed earlier this is independent of the choice of gauge 
on each site. The properties of the system in this temper- 
ature region are clearly very different from Fermi liquid 
behaviour as Anderson has stressed and these two energy 
scales should correspond to his ‘holon’ and ‘spinon’ en- 
ergy scales respectively. The thermopower should obey 
the Heikes formula [24] and we can expect only a low 
mobility of the holes. However, a more detailed study of 
this regime is required. 


VI. DISCUSSION 


for the ground state at the half-filled band is highly re- 
dundant, due to a local SU( 2) invariance at exactly half- 
filling. This redundancy is reflected in an apparent de- 
generacy of the BCS trial wave-function before projec- 
tion. Doping destroys the local SU(2) invariance and 
splits these degenerate states, and we find that the sta- 
ble state upon doping is the d-wave pairing RVB state. 
In this RVB state, electrons are paired even at half- 
filling and it costs an energy of order J to break a pair. 
These pre-existing electron pairs lead to a non-zero su- 
perconductivity amplitude upon doping, and the magni- 
tude of this superconducting amplitude or order param- 
eter is shown to be proportional to the hole concentra- 
tion S when 6 is small. The elementary excitations at 
half-filling are the projected BCS quasi-particle states or 
spinons, with four point zeros on the pseudo-Fermi sur- 
face. 

Our analytic approach can also be applied to ID and 
systems with dimensionality d > 3. We find that lowest 
energy state in ID is the projected Fermi liquid RVB 
without electron pairing, as shown in Appendix 2. Our 
theory predicts no superconductivity in a ID RVB. For 
large d, the energy per bond in the RVB pairing state 
is proportional to 1/d, reduces relative to an AF. So the 
pairing state is particularly favourable in 2D. The precise 
form of the 2D phase diagrams which depends sensitively 
on the relative energies of the AF and d-wave RVB states 
as a function of S is too subtle a question to be settled 
by the approximation we use here. 

There are many questions that require further investi- 
gation such as the exact relationship between the discus- 
sion here in terms of phase coherence among the paired 
electrons and Anderson’s ’holon’ [25] concept or the na- 
ture of the high-temperature phase where this phase co- 
herence is lost but strong magnetic correlations remain 
and presumably do not lead to a Fermi liquid that is the 
usual description of a normal state. 
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We have used a variational method to study a pro- 
jected BCS trial wave-function for the square lattice ef- 
fective Hamiltonian. Using a Gutzwiller approximation 
to treat the effect of the projection operator, we ob- 
tained a renormalised Hamiltonian in which the pro- 
jection operator is replaced by renormalisation factors. 
This approximation is shown to be in good agreement 
with numerical Monte Carlo calculations for such pro- 
jected wave-functions. In this mean-field approxima- 
tion both particle-particle and particle-hole pairing am- 
plitudes must be included. The fermion representation 


VIII. APPENDIX 

A. The projected spin-density-wave state 

In this Appendix, we use the renornralised Hamilto- 
nian approach to study the projected spin-density-wave 
state for effective Hamiltonian (1). That state was pro- 
posed and studied using VMC by Yokoyama and Shiba 
[18]. The generalisation of the Gutzwiller method to the 
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antiferromagnetic states for the hopping process was for- 
mulated by Ogawa and co-workers [26] . 

The projected spin-density- wave state [18], [19] is 

\ip) = Pd IV’o) (Al.l) 

IV'o) = II (UfcC fc,<T + aV k4+Q,a) 1°) ( AL2 ) 

ka 

where k runs over the Fermi sea, Q = 7r/a(l, 1), and 

u k = [(1 + cos0 fc )/2] 1/2 
v k = [(1 - cos6» fc )/2] 1/2 
cos 9 k = lk/(^ 2 AF+ll) 1/2 

A af is a variational parameter, and is given by (18c). 

In a study of the expectation value in the state (Al.l), 
we use the Gutzwiller approximation to replace the pro- 
jection operator by renormalisation factors. In analogy 
to the analysis we discussed in section2, we find that 

1 — n 

1 — 2 n^ni/n 
gs = { 1 - 2n^n l /n) 2 

where n- \ and n± are the spin-up and spin-down electron 
occupation number of state | ipo) in one sublattice respec- 
tively. The renormalisation factors reduce to (4a, 4b) in 
the case n j = , and the form for g t agrees with [26] . 

Within this scheme, we obtain the energy per site at 
half-filling 

w = -2J(a 2 + 6b 2 )/(l + a 2 ) 2 

where 

a = N s 'YAaf/ (A ~ AF + 7fc) 1/2 

k 

b = (SNsr^l/i^AF+ri) 1 ' 2 

k 

The spin-spin correlation (S) • Sj) = | w , and the stag- 
gered magnetisation is 

Ms = rTf — nf = 2a/ (1 + a 2 ) 

with rv, the occupation number in the state \i[f). (Si • Sj) 
and Ms, are plotted in figure Al as functions of A af, 
Figure Al also shows the VMC calculations [18]. The 
case A af = 0 corresponds to the projected Fermi liquid 

state, while Aaf — > oo corresponds to the Neel state. 
The results agree well for large values of A af, but there 
are substantial deviations for small A af- 
The kinetic energy per hole in our theory is 

T = —16tb/ (1 + a 2 ) 

For the optimal value of Aaf{~ 0.9), T = — 2.16f, sub- 
stantially higher than that in the d-wave pairing state. 
Note that this value is also higher than that found in the 
VMC calculation [18], where the optimal AAF is found 
to be much smaller. 



FIG. 9: Spin-spin correlation and staggered magnetisation 
Ms as functions of A af in the projected spin-density-wave 
state. The full curves are the results of the renormalised 
Hamiltonian approach, and the broken curves are the VMC 
results (extrapolated to the infinite systems) by Yokoyama 
and Shiba [18], [19]. 


B. RVB in a ID system 


The renormalised Hamiltonian approach can be 
straightforwardly applied to the model (1) in ID. Using 
the projected BCS wave-function (2), and the same tech- 
nique for 2D, we have found that the RVB ground state 
at the half- filling is described by an equation between C , x , 
and A x (defined in (12)-(13)): 


Cl + 



= C 2 


(A2.1) 


with 

Ci = (2N s )~ 1 '^2 l C0S (MI = 2/7T 2 

k 


(A2.1) is parallel to (17) in 2D. Similar to the 2D case, 
different parameters in (A2.1) are related to each other 
under the SU(2) gauge transformation, and correspond 
to the same physical state. This state is described by the 
projected Fermi liquid state, where A x = 0, C, x = C\. 
Earlier Monte Carlo calculations [27] and more recent ex- 
act calculation [28] with this wave-function have shown 
that the energy of this state is extremely close to the 
exact solution [29]. The energy per site in our analytic 
mean- field approach is — 6/ 7r 2 . This value deviates by 
about 37% from the true result [27], [28]. This quanti- 
tative discrepancy however is not surprising, because the 
Gutzwiller approximation is poor in ID. 

In parallel to the discussion in section4, we can study 
the system with some holes. We found that the stable 
lowest-energy state corresponds to A x = 0. Introducing 
the finite value of A x , the system loses kinetic energy. 
Therefore we expect there is no electron pairing and no- 
superconductivity in this ID RVB. 
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